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f^*) I Abstract. Discontinuities in relativistic hydrodynamics - shock waves and 

freeze-out shocks - are considered across both space-like and time-like hyper- 
surfaces. We analyze the peculiar features of the freeze-out discontinuities and 

their connection to the shock-wave phenomena. 
-*— » ■ 

"o . 

1. Hydrodynamical Discontinuities 

Ph ■ 

Relativistic hydrodynamical model M has been widely discussed in recent years 
within their connection to high energy nucleus- nucleus (A+A) collisions (see, for 
example, |2j]). The system evolution in relativistic hydrodynamics is governed by 
the energy-momentum tensor 

T" v = {e+p)u^u v - pg^ (1) 

and conserved charge currents. The baryonic current, 

f = n u" , (2) 

plays the main role in the application to A+A collisions. 

The hydrodynamical description includes the local thermodynamical variables 
(energy density e, pressure p, baryonic density n) and the collective four- velocity 
u^ = (1 — v 2 ) _1 / 2 (l, v). The continuous flows are the solutions of the hydrodynam- 
ical equations 

d^ v = , d„f = , (3) 
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with specified initial conditions. 

Equations (0) are the differential form of the energy-momentum and baryonic 
number conservation laws. Along with these continuous flows, the conservation laws 
can also be realized in the form of discontinuous hydrodynamical flows which are 
called shock waves and satisfy the following equations: 

T^K=T^A V , no<A M = nu"A^ , (4) 

where A M is the unit 4-vector normal to the discontinuity hypersurface. In Eq. (^) 
the zero subscript corresponds to the initial state ahead of the shock front and the 
quantities without an index are the final state values behind it. To complete the 
system of the hydrodynamical equations (J3|) or (|j) the fluid equation of state (EoS) 
has to be added as an input p = p(e, n). 

Besides the fluid EoS and the initial conditions for the fluid evolution the 
third crucial ingredient of the hydrodynamical model of A+A collisions is the so- 
called freeze-out (FO) procedure, i.e. the prescription for the calculations of the 
final hadron observables. Particles which leave the system and reach the detectors 
are considered via FO scheme, where the frozen-out particles are formed on a 3- 
dimensional hypersurface in space-time. The FO description has a straight influence 
on particle momentum spectra and therefore on all hadron observables. A gener- 
alization of the well known Cooper-Frye (CF) formula p| to the case of time-like 
(t.l.) FO hypersurface was suggested in Rcf. B. The new formula does not contain 
negative particle number contributions on t.l. FO hypersurface appeared in the 
CF procedure from those particles which cannot leave the fluid during its hydro- 
dynamical expansion. The FO procedure of Ref. M has been further developed in 
a series of publications [5-10]. The particle emission from the t.l. parts of the FO 
hypersurface looks as a 'discontinuity' in the hydrodynamic motion. We call it the 
FO shock. 

In the present paper we analyze the discontinuities in relativistic hydrodynamics 
- normal shock waves and the FO shocks - across both space-like (s.l.) and t.l. 
hypersurfaces. 

2. Relativistic Shock Waves 

The relativistic shock waves are defined by equations (0) . The important constraint 
on the shock transitions (m is the requirement of non-decreasing entropy (thermo- 
dynamical stability condition): 

su^A^ > SoU^Af, , (5) 

where s is the entropy density. 

We consider one-dimensional hydrodynamical motion in what follows. In its 
usual sense the theory of the shock waves corresponds to the discontinuities across 
the t.l. surface, i.e. the normal vector A^ is a s.l. one. It means that the shock-front 
velocity is smaller than 1. In this case one can always choose the Lorentz frame 
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where the shock front is at rest. Then the hypersurface of shock discontinuity is 
x a h — const and A^ = (0, 1). The shock equations (0) become 

T oi = T oi ; T n = T u noU i = nu i (6) 

From Eq. (H) one obtains 



v 2 _ JP - Po)( £ + go) ^2 _ (P-Po)(£q+P) / 7 s 

° (e-e )(e +p) (e-e )(e+p ) 

Substituting (|7|) into the last equation in (|6|) we obtain the well known Taub adiabate 
(TA) equation |IJ 

n 2 X 2 - n 2 Q X 2 - (p - Po )(Jf + X ) = , (8) 

in which X = (s + p)/n 2 , and TA therefore contains only the thermodynamical 
variables. 

The point (e ,p ,n ) is called the center of the TA. The mechanical stability of 
the shock transition from the state (e ,p ,n ) to (e,p,n) was studied in Refs. 113, 
^3| . Thermodynamical stability (||) follows from the mechanical stability and the 
inverse statement is not in general true |12,n3[. The consequences of the mechanical 
stability are also the well-known inequalities for the speed of sound and the flow 
velocities at both sides of the shock front in its rest frame 

c so < v , c s > v . (9) 

In the thermodynamically normal media the compression shocks are stable 
whereas the rarefaction shocks become stable in the thermodynamically anoma- 
lous media. The shock-wave stability in the case of the phase transitions between 
hadron matter and the quark-gluon plasma are discussed in Refs. |14|,|l5|. 

Let us consider now the discontinuities on a hypersurface with a t.l. normal 
vector A M (t.l. shocks). This new possibility was suggested by Csernai in Ref. [|16[ . 
In this case one can always choose another convenient Lorentz frame ('simultaneous 
system') where the hypersurface of the discontinuity is t 8 h — const and A„ = (1, 0). 
Equations (|j) become then 

(10) 



rpOQ _ rpQQ rplQ _ rplQ 


n Q u® — nu . 


From Eq. (|l(]) we find 




2 (e-e )(e +p) _ 2 
V ° (p- Po )(e+p ) ' ' 


_ (e - e )(e + p ) 
(p-p )(e +p) 



(11) 

where we use " ~ " sign to distinguish the t.l. shock case (O) from the standard s.l. 
shocks (0). Substituting (^lj) into the last equation in (|lO|) one finds the equation 
for t.l. shocks which is identical to the TA of Eq. flq). We stress, however, that 
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the intermediate steps are quite different. The two solutions, Eqs. (O) and ([?]), are 
connected to each other by simple relations, 



£ 



(12) 



between the velocities of s.l. and t.l. shocks. These relations show that only one 
type of transitions can be realized for the given initial and final states. Physical 
regions [0, 1) for v^, v 2 (0) and for v 2 , v 2 ( pd| ) in the (e-p)-plane are shown in Fig. 1. 

When the initial and final states are thermodynamically equilibrated the TA 
passes through the point (e 0l p ) and lies as a whole in the regions I and IV in 
Fig. 1 p3| , i.e. only compression and/or rarefaction s.l. shocks (with s.l. normal 
vector A^) are permitted. 

The only way to make the t.l. shocks to be possible is to allow the metastable 
initial and/or final states. Then new possibilities of t.l. shock transitions (fLO 11) to 
regions III and VI in Fig. 1 would be realized (see, e.g., the t.l. shock hadronization 
of the supecooled quark-gluon plasma in Ref. Il7|] ) . 
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Fig. 1. Possible final states in the e — p plane for shock transitions from the 
initial state (e ,p ). I and IV are the physical regions for s.l. shocks, III and VI for 
t.l. shocks. II and V are unphysical regions for both types of shocks. Note, that 
only states with p < e are possible for any physical equation of state in relativistic 
theory. 
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3. Freeze Out Shocks 

In the hydrodynamical model of A+A collisions the fluid expansion has to be ended 
by the correct transformation of the fluid into free streaming hadrons. The most 
important requirement for this FO procedure is to satisfy energy-momentum and 
charge conservations during the fluid transition into free final particles. It is assumed 
that there is a narrow space-time region where the mean free path of the fluid 
constituents increases rapidly and becomes comparable with the characteristic size 
of the system. The local thermal equilibrium is supposedly maintained in this 
intermediate region. In practice, one considers a "zero width" approximation and 
introduces a FO surface, so that particle distributions remain frozen-out from there 
on. 

In his original paper Landau [Q defined the FO hypersurface by the condition 
T(t, x) = T* , where t and x are the time and the space coordinate respectively, 
and T* is the fluid FO temperature chosen to be approximately equal to the pion 
mass. The FO procedure, first introduced by Milekhin ]19|j , was improved further 
by Cooper and Frye |3(] and this method has been used ever since. Final momentum 
spectra for i-th type hadrons are expressed by the formula J3[ : 

d 3 N t di f f k v u v -^ \ 

k -w = m 4 d ^ k * \—f^) ■ (13) 

where da m u is the external normal 4- vector to the FO hypersurface £/. It equals 
to Str{—dx, dt) in 1+1 dimension, where St r is the (constant) transverse area of the 
system, fa denotes the local thermal Bose or Fermi distributions, di is the degener- 
acy factor for particle i with the chemical potential /ij. The particle spectrum (O) 
is known as the CF distribution function. 

The initial conditions of the hydrodynamical motion are given on the initial (s.l.) 
hypersurface Y, in . The final hypersurface £/ should be closed to T, in and in general 
£/ consists of both s.l. and t.l. parts. The CF formula (|l3|) leads automatically to 
the energy-momentum and charge number conservations, if the ideal gas fluid EoS 
at the FO hypersurface is assumed. 

The CF formula ( |l3| ) still does not provide a complete solution of the FO 
problem. The FO surface consists of t.l. parts too. Eq. ( ]l^ ) can not however be 
used for a t.l. FO surface (i.e., s.l. normal vector da^): free final particles "return" 
to the fluid if da^k^ < 0, and this causes unphysical negative contributions to the 
number of final particles. The modified FO procedure and new formula for the final 
particle spectra emitted from the t.l. FO hypersurface was proposed in Ref. Q: 

,d 3 N t di f , ,„, , fk x u x Q _fif 



k °^r = ?2# L d ^ * -^ I ^*") • < 14 ) 



Eq. ( |l4| ) looks like CF formula (|13|), but without negative particle numbers that 
appear in (|]j|) for t.l. FO surfaces. These negative contributions are cut-off by the 
6>-function in Eq. (Oh. We'll call Eq. (ph the cut-off (CO) distribution in what 
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follows. An inclusion of the CO FO into the self-consistent hydrodynamical scheme 
is considered in Ref. |nj . The distribution function for the free particles in Eq. ( |l4| ) 
contains the new parameters T g ,fif,u" We'll briefly call the free particle state as 
the 'gas' in order to distinguish it from the ordinary fluid. Note, however, that 
final particle system differs from the normal fluid and gas. It has the non-thermal 
CO distribution function which is frozen-out as all further particle rescatterings are 
'forbidden' in the post-FO gas state. 

The presence of the ^-function in the right hand side of Eq. (p4h leads to the 
discontinuity between the fluid Tj 1 \x\ 1 uK and the gas T g , nf,u g variables across the 
t.l. FO hypersurface to satisfy the energy-momentum and charge conservations. 
We call this discontinuity the FO shock. 

Let us consider this FO shock on t.l. hypersurface in more detail. To obtain 
the analytical solution of the problem we restrict ourself to the case of zero baryonic 
number n — and consider massless particles with the Jiittner distribution function 

4£) _„(_*£). („, 

The energy-momentum conservation between the fluid and free particles along the 
t.l. hypersurface acquires the form: 



" f (2tt) 3 J k° \ T f 

= *>JT - (gfs / ^ *"*" * i^f) '(****) ■ (16) 

The energy-momentum tensor T*t v in Eq. (|16|) has the form of Eq. (|lj) with the 
functions e(T/) = Ef and p(Tf) — pf given by 

e(T) = 3p(T) = JLjfVd** exp(-^ = ± T 4 . (17) 

In the rest frame of the FO front, i.e. in the Lorentz frame where da^ becomes 
equal to S't r .(0, dt), the T® 1 and T 11 components of the gas energy-momentum tensor 
in Eq. (Qq) can be rewritten then in the form 

Tf = (e* g +p* g )uy g , Tf = (e* g +p* g ) U y g +p* g . (18) 

It also coincides with that of Eq. (|l|) , but with effective values of e* and p* which 
are found to be equal to 

4-Pi.M, P ; = P (r 8 ) " + ".y-.» , ( 19 ) 

where v g is the velocity parameter of the gas in the FO shock rest frame. The 
functions e and p in the right hand side of Eqs. (19h are given by Eq. (fij 
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Due to the same formal structure of r°\T g n @ and Tf, Tj 1 (Q) one obtains 
the solution of Eq. ( |l6| ) 

2 (p/-Pp(e;+P/) 2 (p/-p;)(g/+y;) , m 

which is similar to Eq. (pp. Note, however, that in contrast to Eq. (m) the values of 
e* and p* in Eq. (EG) are not just the thermodynamical quantities, but depend also 
on v g . 

It should be emphasized that Eq. fll8| ) can also be introduced for the distribu- 
tion functions of massive and charged particles [|10| . In the latter case the energy- 
momentum conservation leads to the familiar expressions ( |20|) for the effective en- 
ergy density and pressure, and the charge conservation leads to the TA equation 
(||) for the effective charge density. 

We fix the FO hypersurface £/ by the condition T g — const. The analytical 
solution of Eq. ( |20| ) can be presented then in the form 

9v} -8v f + 3 _ (T f \ A _ ^v}-2v } + lf(v s + l) 



Sv} + 1 R ~ [tJ (3v f -l)(3v} + iy ■ (21) 

It is instructive to compare these results with shock-wave solution tm) for the same 
ideal gas EoS (|l7|). We fix the temperature T g of the final state in s.l. shock wave 
(pi) and consider v g and Tf dependence on vf. 

1 n (T f V 3(1 -vj) 



Figs. 2 and 3 show the dependences of v g and Tf on Vf for the fixed temperature 
T Q in the final state, both for the FO shock ( |2l| ) and for the normal shock wave 
(f22|). The kinematical restrictions on the gas velocity give the same value of the 
minimal fluid velocity in both shock transitions, (vf) m in — pf/ef. It equals to 1/3 
for the considered ideal gas EoS ( J17| ) of the fluid. Figs. 2 and 3 indicate that at 
low values of the fluid velocity, Vf < l/\/3, the behavior of Vf and Tf/T g for the 
FO shock (Efl) and for the normal shock wave (E3) is quite similar. The value of 
Vf = l/v3 corresponds to the speed of sound, c s — l/\/3, in the system with ideal 
gas EoS ( |l7| ) . According to the requirements given by Eq. (|S|) the shock transitions 
at Vf < c s are mechanically unstable. Mechanically stable solutions at Vf > c s for 
normal shock waves and FO shocks have qualitatively different behavior. For the 
stable normal shock wave one has v g < Vf and Tf < T g (see Figs. 2, 3), i.e. only 
compression normal shock wave transitions '/ '—>'<? ' would be stable. For the FO 
shock transitions we find a completely different behavior, v g > Vf and Tf > T g , 
illustrated in Figs. 2, 3. In contrast to the normal shock waves, only the rarefaction 
FO shock transitions '/ '— >'<? ' are stable. This result is in agreement with an 
intuitive physical picture of the FO as a rarefaction process. 
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Fig. 2. Gas velocity v g as a function of the fluid velocity Vf in the rest frame 
of the shock front. The dashed lines represent mechanically unstable transitions, 
whereas the solid lines show the mechanically stable FO shocks and normal shock 
waves. Thin lines correspond to a normal shock wave (Eq. (f22[)) and the thick 
ones correspond to the freeze-out shock (Eq. (|2l|)). The circle has the coordinates 
(l/v3; 3 — 4/v3) and the square has the coordinates (1/V3; 1/vo). 



The entropy flux of the fluid is given by 



S^da-ft 



SfU^da^ , 



(23) 



where Sf = (e/ +pf)/Tf is the fluid entropy density. Similar expression is valid for 
the final 'gas' state in the case of normal shock waves. The entropy production in 
the normal shock waves is given by the formula 



s^da^ 
s^da^ 



1 



3 3 / 4 w/ 



9vj 



1 



7 



1/4 



(24) 



and is shown in Fig. 4. The entropy flux of the gas with the cut-off distribution 
function is given by pOl : 



s^da^ = 



day 

(27T)' 



IT** 

k° 






1 



hi 



Ku v g 



The entropy production in the FO shocks is calculated then as 



Sjda^ 



1 

2^7 



(3^ + 1)2(3^-1) 



1/4 



0{da x k x ) . (25) 



(26) 
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The maximal entropy production in the FO shock (Eq) corresponds to Vf = c s = 
1/V3 and can be considered as an analog of the Chapman- Jouguet point (see 
Fig. 4). 




0.5 - 



Fig. 3. Ratio of the temperatures on the both sides of the shock front, Tf/T g , as 
a function of the fluid velocity Vf in the rest frame of the shock front. The legend 
corresponds to the Fig. 2. The circle has the coordinates [l/-\/3; (8/(3y / 3)) 1//4 ] and 
the square has the coordinates (l/\/3; !)• 



4. Conclusions 



In the present paper the discontinuities in relativistic hydrodynamics have been 
considered. The particle emission from t.l. parts of the FO hypersurface looks as a 
discontinuity in the hydrodynamical motion (FO shocks). The connection of these 
FO shocks to the normal shock waves in the relativistic hydrodynamics has been 
analyzed. The above consideration can be used in the hydrodynamical approach to 
the relativistic A+A collisions. It should be also important in the models of A+A 
collisions which combine hydrodynamics for the early stage of the reaction with a 
microscopic hadron description of the later stages (see, for example, Ref. pQ] ). 
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Fig. 4. Ratio of the entropies on the both sides of the shock front, as a function 
of the fluid velocity Vf in the rest frame of the shock front. The legend corre- 
sponds to the Fig. 2. The circle corresponds to the maximal entropy (analog 
of the Chapman- Jouguet point) in the freeze-out shock and it has the coordinates 
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and the square has the coordinates (l/v3; 1). 
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